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Experimental evidence has recently put the validity of the Porter-Thomas distribution (PTD) for 
partial neutron widths into question. We identify two terms in the effective Hamiltonian that violate 
orthogonal invariance (the basis for the PTD). Both are due to the coupling to the decay channels. 
We show that realistic estimates for the coupling to the neutron channel and for non-statistical 
gamma decays yield significant modifications of the PTD. 
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Introduction. Recent experimental results on the dis¬ 
tribution of neutron resonance widths have cast serious 
doubt on the validity of random-matrix theory (RMT) in 
nuclei. RMT predicts that the reduced neutron widths 
follow a Porter-Thomas distribution (PTD) [i| (a dis¬ 
tribution with a single degree of freedom). That predic¬ 
tion assumes non-overlapping resonances with a single 
open channel (the neutron channel). Deviations from the 
PTD for several open channels are analyzed in Refs. Si- 
Agreement with the PTD could be excluded with very 
high probability both for the Pt isotopes Q and for the 
Nuclear Data Ensemble Q. Numerous theoretical at¬ 
tempts @-01 to account for the results of Refs. Hi have 
not definitively resolved the issue. The validity of RMT is 
of central importance for the statistical theory of nuclear 
reactions @| that is widely used in nuclear cross-section 
calculations. 

In the present paper we address two dynamical effects 
that modify the PTD and that apparently have not been 
taken into account so far in the theoretical literature or in 
the analysis of neutron resonance data. The two effects 
are the Thomas-Ehrman shift known from the study of 
light nuclei 0 , and non-statistical effects in the gamma 
decay of the neutron resonances [lH . We show that both 
may cause significant deviations of the distribution of 
neutron resonance widths from the PTD. 

The PTD follows from the orthogonal invariance of the 
Gaussian Orthogonal Ensemble of random matrices (the 
GOE). The GOE distribution function is proportional to 

r A I ^ ^ 

dOexpj . (1) 

^ p /i<IZ (7 


dent. For A —>■ oo, the projections of the eigenfunctions 
onto an arbitrary vector in Hilbert space possess a Gaus¬ 
sian distribution, and the reduced widths, therefore, have 
a PTD. The reported disagreement of the distribution 
of reduced neutron widths with the PTD directly chal¬ 
lenges the postulated orthogonal invariance of the GOE. 
Conversely, dynamical effects that violate orthogonal in¬ 
variance will cause deviations from the PTD. We show 
that both, the Thomas-Ehrman shift and non-statistical 
gamma decays, have that property and, thus, qualify as 
causes of the observed disagreement. 

Violation of Orthogonal Invariance. We consider s- 
wave neutron scattering on a spin zero target nucleus. We 
take account of the single open neutron channel and of 
the large number of gamma decay channels. The effective 
Hamiltonian is @] 

( 2 ) 
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Here E is the neutron energy. We have replaced the 
actual Hamiltonian by the GOE Hamiltonian iy^OE g^g 
defined in Eq. O- The real matrix elements {E) with 
/i = 1, 2,..., A couple the s-wave neutron channel to the 
space of A resonance states and carry the same energy 
dependence as do the neutron partial width ampli¬ 
tudes. The shift matrix E accounts for that energy de¬ 
pendence. It is the analog of the Thomas-Ehrman shift, 
with elements 


^OO 

F^,{E) = Pv / dE 
Jo 


„ W ^iE')W,iE') 
E-E' 


( 3 ) 


Here A with A —oo is the dimension of the GOE ma¬ 
trices. The parameter A defines the width 4A of the GOE 
spectrum. The mean level spacing at the center of the 
GOE spectrum is d = ttA/A. The E^ are the GOE eigen¬ 
values, and dO is the Haar measure of the orthogonal 
group in A dimensions. It encompasses the GOE eigen¬ 
functions. Factorization of the distribution implies that 
eigenvalues and eigenfunctions are statistically indepen¬ 


where Pv indicates the principal-value integral. At ener¬ 
gies far above threshold the matrix E is often neglected 
because then contributions to the integral from energies 
E' < E and E' > E tend to cancel. Such cancellation 
cannot occur at neutron threshold A = 0 and E^v{E) 
may, thus, not be negligible. We neglect contributions 
similar to A from closed channels and take A(A) as a 
paradigmatic example. The matrix elements play 
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the same role for the gamma channels as do the for 
the neutron channel except for a different energy depen¬ 
dence of the , resulting in a negligible contribution 
to the Thomas-Ehrman shift. 

None of the terms added to in Eq. © is invari¬ 

ant under orthogonal transformations. Addressing the 
regime of isolated resonances we confine our attention to 
the matrix F and to the coupling to the gamma decay 
channels. Thus, the elements of the width matrix W^Wi, 
only serve to define the neutron decay widths and are 
otherwise negligible. Beyond that regime, the width ma¬ 
trix does cause deviations from the PTD [3,0- 

Let us first disregard gamma decay channels and con¬ 
centrate on the role of F. Writing W^{E) = 
we observe that the matrix has a single nonzero 

eigenvalue The associated eigenvector defines 

the superradiant state 13, 3| labeled /r = 1. The trans¬ 
formation to the eigenvector basis leaves the ensemble of 
GOE matrices unchanged and yields 


h: 


eff 


fiiy 




F-inJ2w^{E) 


(4) 


where F = Y.^ F^^. 

The transformation does not diagonalize the matrix F 
exactly because the integral defining F receives contribu¬ 
tions also from higher energies where the approximation 
W^{E) K. does not apply. Near neutron thresh¬ 

old such contributions are relatively small, however, and 
Eq. ([4]) should be a good approximation. 

We compare F (Eq. (|3])) with the diagonal element 
Sufficiently far above neutron threshold (where 
the s-wave penetration factor is ~ 1) the matrix ele¬ 
ments Wfj, obey [3| = Ndx/n'^ where x is re¬ 

lated to the average of the s-wave scattering function via 
S' = (1 — x)/{l + x) and is, thus, of order unity. Hence 
~ dN/TT^ = A/tt. For N ^ 1 that is much larger 
than the root-mean-square value \^J2/N of the diagonal 
element For A —>■ 0 the s-wave penetration fac¬ 
tor reduces The reduction does not affect the 

principal-value integral, however, which actually attains 
its maximum value at A = 0. Thus, we expect F Ri A/tt. 
The effect of F can be amplified beyond our estimate by 
a single-particle resonance near neutron threshold. Such 
is the case in the Pt isotopes for the 4s-state of the shell 
model [1]. 

Next, we turn to the gamma channels in Eq. ([S]). Ac¬ 
cording to the statistical model, the matrix elements 
are Gaussian-distributed random variables. The 
total gamma widths of the neutron resonances are, there¬ 
fore, expected to have a distribution with a large num¬ 
ber of degrees of freedom, and the partial widths P^^^ of 
the neutron resonances in a given target nucleus are ex¬ 
pected to have all nearly the same value. The effective 
Hamiltonian is, thus, expected to be approximately given 


by - iT^W^Wv - 5fji,viT/2, with P indepen¬ 
dent of ji. The term is obviously orthogonally 

invariant. While the contribution of each gamma chan¬ 
nel to r is small the number of such channels is large 
resulting in a value of T that dominates the total neu¬ 
tron resonance widths near neutron threshold. 

A recent analysis ll| of the distribution of total 
gamma decay widths of neutron resonances in ®®Mo con¬ 
tradicts the expectation that these all have the same 
value. For s-wave resonances and positive parity states, 
the distribution is shown in the lower part of Fig. 6 of 
Ref. ll|. The distribution is much wider than predicted 
by the GOE. The result confirms earlier data [15| com¬ 
prising a much smaller number of resonances. It seems 
that at present, the cause for the deviation is not un¬ 
derstood. It is not clear whether it is due to a specific 
property of ®®Mo or whether it is likely to occur in other 
nuclei as well. We opt for the second possibility. We as¬ 
sume that the total gamma decay widths generically pos¬ 
sess large fluctuations, in contradiction to the statistical 
model. We explore the consequences of that hypothesis 
for the distribution of neutron decay widths. 

Typical values of total gamma decay widths for s- 
wave neutron resonances are of order 100 meV, both in 
medium-weight 111 and in heavy [l6| nuclei. The dis¬ 
tributions for the total gamma decay widths shown in 
Ref. [HI start roughly at To = 100 meV and fall off with 
a half width a of roughly 300 meV. To transcribe these 
figures into the effective Hamiltonian of Eq. m we use 
that in medium-weight and heavy nuclei typical average 
resonance spacings d = ttX/N are of order 10 eV. Then 
To Ri ttA/ (lOOA^) and a m ttA/ (30A^). The orthogonal in¬ 
variance of the coupling to the gamma decay channels is 
broken by the spread a. We compare a with the Thomas- 
Ehrman shift function F r; A/tt. We have (l/A^)Tr a r; 
7rA/(30iV) while (l/Af)Tr = F/N = A/(7riV). 

While somewhat smaller, gamma decay breaks orthogo¬ 
nal invariance roughly as strongly as does the Thomas- 
Ehrman shift. 

In summary we have identified two effects that break 
the orthogonal invariance of the GOE, the Thomas- 
Ehrman shift and the non-statistical distribution of 
gamma decay widths. In the remainder of the paper we 
investigate the consequences of both effects for the dis¬ 
tribution of neutron resonance widths. We do so in the 
framework of a schematic model. We write the effective 
Hamiltonian as 


= H^^^+d^iS,iZ. (5) 

We have suppressed the constant gamma decay width T 
because it only causes a uniform shift of all eigenvalues 
[3] and does not affect the eigenfunctions. The constant 

Z = F - {i/2)dT - iTT (6) 
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violates orthogonal invariance and includes the Thomas- 
Ehrman shift; the effect of non-statistical gamma decays 
where we schematically represent the spread of gamma 
decay widths by a single term; and the term —Wp 
which is a reminder that we use a basis where the state 
|1) is the superradiant state. 

Average Level Density. It is easy to see that for N —>■ 
oo the presence of Z in the effective Hamiltonian has 
a negligible influence on the average level density. The 
level density is defined as p{E) = —(l/7r)S(£’+ — 

Here (E+ — is the retarded Green function. This 

expression for p is physically meaningful only if i?®® is 
Hermitian, i.e., if Z is real. We first consider that case 
and put Z = F. The term F in H®® may be consid¬ 
ered as causing a doorway state. Therefore, we treat the 
first line and the first column of H®® differently from 
the rest. We define the orthogonal projection operators 
P = |1)(1| and Q = 1 — P. With > 2 we write 

PiJ®®P = PpOE + z = Fo, (PP®®g)i^ = = Vf,, 

(QP®®g)m = PGOE ^ The el¬ 

ements of the matrix QH^^Q are Gaussian-distributed 
random variables. Moreover, the probability distribution 
of QH^^Q is invariant under orthogonal transformations 
in g-space. Therefore, the matrices QFl’^^Q form a GOE 
of dimension — 1. We denote that ensemble by 
and, suppressing the term J2p i write the total Hamil¬ 
tonian in matrix form, 

rrefF _ f ^0 ^ \ 

The right-hand side of Eq. 0 is identical in form with 
the standard model for a doorway-state (see, for instance, 
Ref. [13 )■ The spreading width of the doorway state is 
rl = 27r(l/A^) Vjf/d. In the present case we have 
Vp = Sp Hy. We note that for TV 1 the sum is 
self-averaging and given by for every member of the 
ensemble ©. The resulting value of the spreading width 
r-l = 2A is comparable with the total width of the GOE 
spectrum, and the doorway state is completely smeared 
over that spectrum. While in the standard doorway- 
state model [ 13 , the magnitude of the coupling matrix 
elements is of order d, here their root-mean square 
values d'/N are fixed by the underlying GOE and very 
large compared to d. Therefore, the doorway state does 
not cause a local enhancement of the level density. We 
note that the model of Eq. 0 is physically meaningful 
only for jZj < 2A. For jZj ^ 2A there exists a distinct 
state outside the GOE spectrum that carries (almost) 
all the coupling to the neutron channel. That case does 
not seem to model the scattering of slow neutrons in a 
meaningful way. 

As a corollary we mention that for imaginary Z and 
jZj Ri 2A we deal with a superradiant state that causes a 
pole of the Green function in the complex energy plane. 
The distance of that pole from the real axis is comparable 


with the width of the GOE spectrum. Therefore, the 
superradiant state conveys only a small part of its large 
neutron width to the remaining GOE eigenstates. It does 
not seem meaningful to consider imaginary values of Z 
in excess of 2A. 

Numerical Results. With (j)p the normalized eigenfunc¬ 
tions of the effective Hamiltonian ([5]) with J2p ^p 0, 
the partial neutron decay widths are proportional to 
A^|(l|^p)p. The factor N is introduced so that the av¬ 
erage width equals unity. We have calculated the effect 
of Z on the distribution of partial widths x = A^|(l|^p)p 
perturbatively, both for \Z\ <C d and for \Z\ ^ A. The 
results show that Z does influence the PTD. Both limits 
are unrealistic, however, and serve only as a check for the 
numerical work. Therefore, we do not give our analytical 
results here. 

We first present results for real Z. We use the dimen¬ 
sionless parameter k = ZjX and consider values that are 
physically realistic but lie outside the range of validity 
of the perturbative approach. As Z increases the PTD 
Plx) = (\j(2\/2'Kx')') exp{—a:/2} is deformed. The re¬ 
sulting probability distribution V{x) of the partial widths 
is shown in a plot of \P(x)/P{x)\ — I versus x for several 
values of k in Fig. [T] The term Z'|l)(l| in iJ®® leads to 
a segregation of states. The states in one group become 
broader and those in the other group become more nar¬ 
row. Therefore, the modified distribution has a longer 
tail and is more strongly peaked at x = 0 than the PTD, 
while in the middle around x = I the distribution is sup¬ 
pressed. In the limit of very large k one state becomes 
collective (i.e., carries almost all the decay strength). The 
distribution of widths of the remaining states returns to 
the PTD but with a much reduced average width. That 
is consistent with our perturbative results for large Z. 



FIG. 1: (Color online) Relative difference \P{x)/P{x)\ — 1 
as a function of x for several values of k as indicated in the 
figure and for N = 1000. The dash-dot lines are fits using 
expression ®. 

We have fitted these curves using the parametrization 

V{x) = [I + A(1 — x) -I- B{x^ —Qx + 3)]P(x) (8) 

suggested by the perturbative result for small Z. The 
distribution V{x) in Eq. (I5|) is normalized to unity for 
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all values of A and B. Terms linear (quadratic) in x 
embody a change in average width (a quenching of the 
PTD, respectively). For each of the curves representing 
the data we also show a fit using Eq. dS]) as a dash-dotted 
line. In the interval 0.1 < k < 1 that is most relevant for 
applications, the deviations from the PTD are described 
perfectly by expression ([8]). 

The scaled fit coefficients A/k? and B /are shown 
as functions of k in Fig. [2] for different matrix dimension 
N. For K < 0.4 we encounter numerical instabilities. 
However, in that regime the changes of the PTD are too 
small to be of practical interest. For k > 0.4 the results 
for different matrix dimension N exhibit consistently a 
linear dependence on k approximately given by 

A/k^ = -0.035 ± 0.010-b (0.16 ±0.01 )k , 

B/k^ = 0.146 ± 0.002-b (0.099 ± 0.003 )k . (9) 



0 0.2 0.4 0.6 0.8 1 


K 

FIG. 2: (Color online) Scaled fit coefficients A/t^ and B/k^ 
with A and B defined in Eq. are plotted as functions of 
K for GOE ensembles of different dimension N as indicated 
in the figure. For N = 1000 the error in the fit is shown by 
error bars. The errors for the other curves are similar. The 
double-dot-dash black line shows the linear fits ® of 
and B /as functions of ac. 

As a single quantitative measure of deviations from the 
PTD we use the coefficient of L variation (also known as 
the Gini coefficient) defined as 

T = — J dx J dx' \x — x'\'P{x)'P{x'). (10) 

The advantage of using r over the traditional coefficient 
of variation is that r is nearly insensitive to the existence 
of a small fraction of highly collective states that may 
comprise the average value of a; Q. The coefficient r 
ranges between zero and unity. For a distribution 


with iz = 1 and v = 2 we have r = 2/7r and r = 1/2, 
respectively; r decreases with increasing v. A value of 
T > 2/7r Ri 0.64 corresponds to a distribution that is more 
strongly peaked at small x than the PTD and effectively 
has < 1. For a distribution of the form m we find 
T = (2/7r)(l -2A-A^ + 2B-2AB - W^). For k = 0.8, 
for example, Eq. dH) yields r = 0.70, in agreement with 
numerical tests. That value of r has a x^ distribution 
with v = 0.72. 

The results for purely imaginary Z are qualitatively 
similar to those for real Z. Compared to the PTD, the 
distribution is increased for small and large x and is de¬ 
pressed for a; Ri 1. That is shown for ZjX = —0.8f by 
the black solid line in Fig. [TJ For imaginary Z we do not 
present a fit because we have not succeeded in finding 
similarly good fit formulas as in Eqs. m for real Z. Con¬ 
clusions. The PTD follows from the orthogonal invari¬ 
ance of the GOE. We have identified two causes for vio¬ 
lation of that invariance: The Thomas-Ehrman shift and 
non-statistical gamma decays. We have shown that rea¬ 
sonable estimates for both cause significant deviations of 
the distribution of neutron decay widths from the PTD. 

Invariance breaking by the Thomas-Ehrman shift is 
due to the coupling to the neutron channel. Such cou¬ 
pling is immanent in the theory and does not invalidate 
the GOE. The shift is strongest at neutron threshold and 
is expected to be particularly pronounced when the s- 
wave strength function is maximal. In contradistinction, 
it may be argued that the existence of non-statistical 
gamma decays represents a genuine violation of GOE as¬ 
sumptions. It is conceivable that such decays are due 
to transitions to low-lying states where random-matrix 
theory does not apply. Hopefully time will tell. 

In all cases studied, invariance breaking results in a 
depletion of the probability distribution for the partial 
widths X near x = 1, compensated by an increase for 
small and large values of x. To estimate the effect of 
such breaking in an individual nucleus, the quantities F 
and dr in Eq. ([51) must be estimated. For F that should 
be possible using the neutron strength function. For dF 
(a measure of the spread of total gamma decay widths) 
the simultaneous analysis of the distribution of partial 
neutron widths and of total gamma decay widths for the 
measured chain of neutron resonances is required. 
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